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In these proceedings, we shall first recall the evolution equations arising when 
increasing the rapidity Y = log(s) within the perturbative QCD regime. We shall 
then summarise the main properties on their asymptotic solutions and discuss the 
physical picture emerging from our analysis. 



We start with the equations describing the evolution of scattering am- 
plitudes towards high-energy. For simplicity, we shall work in the large- N c 
limit for which we can use the dipole picture. Thus, we consider the scat- 
tering amplitude (T^ yi ;...;x k y h ^ between a projectile made of fast-moving 
colourless qq dipoles of transverse coordinates (x^y,) and a generic target. 

The evolution equation for (T xy ) can be obtained from dipole splitting 
in the projectile (the large- N c version of a gluon emission) which directly 



leads to 

Q /"ja. ( x -y) 2 



dv <T xy ) = £ J <Pz (x _^ )2( ^_ y)2 [(T xz ) + (T zy > - (T xy ) - (r^ zy ) 

(1) 

The three linear terms in this equation correspond to the well-known BFKL 
equation. This has solutions growing like exp(wpF) (with up = 41og(2)a) 
and thus violating the unitarity bound T xy < 1. The last term, taking into 
account the scattering on both dipoles, is important when the scattering 
becomes of order 1 and restores unitarity. To deal with (1), one usually 
adopts in the mean- field approximation ^T XZ : Zy ^ = (T xz ) (T zy ) which gives 
a closed equation known as the Balitsky-Kovchegov (BK) equation. 

In general, one has to consider an infinite hierarchy of equation for 
/y( fe )V From dipole splitting in the projectile, one gets a BFKL-type 
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contribution proportional to (T^) and saturation corrections going like 
(T( fe+1 )). The resulting infinite set of equations is known as the (large- N c ) 
Balitsky hierarchy. Nevertheless, it has been shown recently that addi- 
tional contributions corresponding to gluon-number fluctuations have to be 
included. Those can be computed through dipole splitting in the target and 
give an additional contribution proportional to ^T^ -1 ') to the evolution 
of (T^ fe '). For the second equation in the hierarchy, this is 

f Muvz^o(xiyi|uz)A(x 2 y 2 |zv)V2 V2 (T uv ) + (1 <-» 2), 

J uvz 

where M. is the dipole splitting kernel and Aq the dipolc-dipole scattering 
amplitude. The new term becomes comparable to the BFKL contribution 
when T ~ a 2 i.e. in the dilute regime as expected. Also, by combining 
pomeron splittings and mergings we can build pomeron loops. Note finally 
that, even if one has to pay the price of an additional factor a 2 , once we have 
a pomeron splitting the amplitude grows like two BFKL pomcrons. This 
becomes comparable to the one-pomeron-exchange for Y > ^log(l/a|). 
We shall discuss in more details the consequences of the fluctuation term 
later on. 

Before considering the solutions of those equations, let us quote that 
the full set of equations can be seen 1 as a reaction-diffusion process in 
which one dipole (x, y) can split into two dipoles (x, z) and (z, y) through 
usual BFKL splitting and two dipoles (xi,yi), (x 2 ,y2) can merge into one 
dipole (u, v) through a vertex directly extracted from (2). However, the 
merging vertex is not positive definite, which proves that this dipole model 
is only effective and fluctuations really involve the quantum gluonic degrees 
of freedom. 

Let us now discuss the asymptotic solutions of those equations. We shall 
neglect impact parameter and work in momentum space k. We start with 
the BK equation which has been shown to lie in the same universality class 
as the Fisher-Kolmogorov-Petrovsky-Piscounov (F-KPP) equation known 
since seventy years in statistical physics. This equation admits travelling 
waves as asymptotic solutions. In terms of QCD variables, this translates 
into the geometric scaling property, stating that T(k, Y) is a function of the 
ration k/Q s (Y) only, where Q S (Y) is the saturation scale which provides a 
natural infrared regulator. The analogy with the F-KPP equation predicts 
that T(k,Y) ~ [k 2 / Q 2 S (Y)]-^ for k > Q s , with Q 2 S {Y) = k 2 exp(av c Y). 
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7 C and v c are pure numbers determined from the BFKL kernel only. It is 
important to notice that the property of geometric scaling is a consequence 
of saturation which extends (as VY) far above the saturation domain, where 
T< 1. 

If we switch on the fluctuation term, the situation becomes a bit more 
intricate. However, neglecting the impact-parameter dependence, one can 
perform a local-noise approximation* 1 which allows to rewrite the infinite 
hierarchy for average amplitudes as a single Langevin equation for event- 
by-event amplitudes which is, formally, the BK equation supplemented by 
a noise term 

d Y T(k, Y) = aK BFKh ®T(k, Y)-aT 2 (k, Y)+a v / 2 K a 2 T(k 1 Y)v(k 1 Y), (3) 

where k is an undetermined constant coming from the local-noise approxi- 
mation and v(k,Y) is a Gaussian white noise satisfying 

(v(k,Y))=0 and Mk, Y)v{k' , Y')) = — S (\og(k 2 /k' 2 )) 6(Y-Y'). 

air 

Again, most of our analytical understanding of the QCD problem comes 
from the fact that (3) lies in the same universality class as the stochastic 
F-KPP equation. The analytical study of that equation in the weak-noise 
limit as well as numerical simulations 3 of (3) show that the effect of the 
fluctuations is double: first, the event-by-event amplitude still displays ge- 
ometric scaling but with an exponent of the saturation scale smaller than 
v c . Then, if one considers a whole set of events, one observes a dispersion of 
their position in \og(k 2 ). This means that the logarithm of the saturation 
scale fluctuates from one event to another, with a dispersion increasing as 
VY, as expected from a random-walk process. A direct consequence of this 
dispersion is that, when we average over all realisations over the noise to 
get the physical amplitude, it leads to violations of geometric scaling. 

To gain more insight on the physical consequences of these results, it is 
sufficient to consider an event-by-event amplitude (r is a reference scale) 

T CT (r, Y) = { 1 9 ~ Ps , p = log ( V 4) and p s = log(r 2 Q% 

which satisfies unitarity constraints, colour transparency and geometric 
scaling, together with a Gaussian probability distribution P(p s ) of aver- 
age (p s ) = log(roOs) an d dispersion a oc Vy. 



a A study of the impact-parameter averaged equations beyond the local-noise approxima- 
tion is under study using a properly defined matching between a mean-field BK evolution 
close to the unitarity limit and random dipole splitting in the dilute tail 2 . 
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The average amplitude is then computed through 

/oo 
d Ps P( Ps )T cv (ps-(Ps)). (4) 
-oo 

At intermediate energies (a < 1), the dispersion of the events is negligible 
and (T((p)) w T cv (p s — (p s )) satisfies the geometric scaling property. How- 
ever, at high energies such that u » 1, dispersion starts to dominate and 
geometric scaling gets violated. A direct computation of (4) shows that, at 
high energy, one has 

™>=<™)=M^) <»'-<*>«' 2 

Thus, within a window increasing linearly with rapidity, it emerges 4 a 
new scaling, called diffusive scaling. Within this window, event dispersion 
becomes larger than the typical decrease length of T ev . We can thus consider 
that, event-by-event, the amplitude is simply a Hcaviside theta function 
<d(p s —p). In other words, at very high energies, the dominant contributions 
to the physical amplitude fully come from fronts which are at saturation i. e. 
by black spots. As an additional consequence, we have (T 2 ) = (T) at high 
energy while the mean field approxiamtion predicts (T 2 ) = (T) 2 . Note that 
all those results can be derived analytically 5 by assuming a strong noise in 
the equation (3). 

Finally, we have to quote that the diffusive scaling property extends (as 
Y) above the average saturation momentum. The fact that we are sensitive 
to saturation, even when the physical amplitude is much smaller than one, 
can have important consequences for LHC physics e.g. for forward jet 
production 6 . 
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